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The magnetic field distribution, the magnetic flux, and the free energy of an 
Abrikosov vortex loop near a flat surface of type-II superconductors are calculated 
in the London approximation. The shape of such a vortex line is a semicircle of arbi- 
trary radius. The interaction of the vortex half-ring and an external homogeneous 
magnetic field applied along the surface is studied. The magnitude of the energy 

■ barrier against the vortex expansion into superconductor is found. The possibilities 
I of formation of an equilibrium vortex line determined by the structure of the applied 

■ magnetic field by creating the expanding vortex loops near the surface of type-II 
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superconductor are discussed. 



1. INTRODUCTION 



X 

I In a mixed state a magnetic field penetrates into type-II superconductors (SC) in 

the form of Abrikosov vortices having one quantum of magnetic flux $o = T^tic/e The 
presence of these vortices in SC and their interaction with inhomogeneities and defects of 
the material (pinning) determine both the magnetic properties of SC and the ability of 
SC to carry a current without a resistance 0. In the vortex core, where the magnetic field 
achieves its maximum, superconductivity is destroyed and the order parameter describing 
the superconducting properties of the material vanishes. The center of the vortex core 
defines the location of the vortex line (VL). The equilibrium form of the vortex lines in 
type-II superconductors reproduce the structure of an applied magnetic field. Thus, the 
homogeneous magnetic field generates a two-dimensional lattice of rectilinear vortex lines 
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directed along the field. The vortex penetration in SC is delayed significantly due to the 
presence of the surface Bean-Livingston (BL) barrier 0, An increase of an applied 
magnetic field Hq results in a decrease of the BL barrier. The barrier is destroyed, 
then the field Ho reaches the thermodynamic critical field Hem- If the applied magnetic 
field Hq < Hem, the penetration of vortices in SC occurs as a result of sufficiently large 
fluctuations. In this case the vortex line produced near the surface forms a semicircle 
(vortex loop) the ends of which touch the superconductor surface P, ^ §]. A minimal 
critical size of the loop exists so that the vortex does not collapse at the surface but 
expands into the SC. 

In this paper we consider the conditions of formation of the vortex loops near a fiat 
surface of a type-II superconductor in the applied magnetic field parallel to the surface. 
In section 2, the magnetic field, the magnetic flux and the free energy of the solitary 
Abrikosov vortex in the form of a semicircle of arbitrary radius are calculated in the 
London approach. In section 3, the interaction of the vortex half-ring and an applied 
homogeneous magnetic field is studied. In this section we calculate the Gibbs free energy 
and determine the energy barrier preventing creation of expanding vortex loops. In section 
4 the processes of formation of equilibrium VLs, defined by the applied magnetic field are 
discussed. The development of a rectilinear VL parallel to the applied magnetic field by 
intersection and cross joining of adjacent vortex loops is considered. 



2. VORTEX HALF RING NEAR THE SURFACE OF A SUPERCONDUCTOR 

Let us consider a type-II superconductor occupying the region z < and bounded by 
the surface S coincident with the coordinate plane (xy). We assume that the vortex line is 
a semicircle of arbitrary radius lying in the plane (yz) parallel to the applied magnetic 
field Ho (fig. IH). Due to the continuity of the magnetic flux, the end points of the VL have 
to be placed at the surface 5* 0. The London approach is appropriate for extreme type-II 
superconductors with a large Ginzburg-Landau parameter k = X/^ ^ 1, where A is the 
magnetic field penetration depth and ^ is the coherence length. In the London model the 
magnetic field distribution H in SC {z < 0) is described by the London equation 0: 
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H + rotrotH = e^(5(r-r„). (1) 

Here is the radius vector defining the VL location in space, and e^, is the unit vector 
directed along the VL. Outside the superconductor {z > 0) there are no currents, and the 
magnetic field H^, produced by the fiux line, is yielded by the Maxwell equations: 

rot H, = 0, div H, = 0. (2) 

The boundary conditions for the surface S can be obtained from the continuity of all 
components of the magnetic field in the plane z = 0: 

H = H,. (3) 

Everywhere here we use dimensionless units and assume that the magnetic fiux $ and the 
coordinate vector r are scaled by the fiux quantum $o and the magnetic field penetration 
depth A, respectively. For these units the scale of the magnetic field is $0/-^^) and the 
current density j is measured by the units of c^o/X^. 

Due to the linearity of equation (|l|) and boundary conditions (|^) we represent the 
magnetic field H in the SC as a superposition: 

H = + Hrf. (4) 

The term H^, = if^e^, describes the structure of the magnetic field of a toroidal Abrikosov 
vortex with a VL consisting of a circle of radius p^. In the cylindrical coordinate system 
(p = \/y'^ + z'^, (f, x) the component H„{p,x) satisfies the equation |]10| : 
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The solution of equation (^) for arbitrary vortex radius p^, was obtained in |10, 11, 12 
by means of the Fourier-Bessel transform, and can be represented as follows: 



H^{p,x) = — / dgexp(«gx)< " " , (6) 
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where /i(C); -^i(C) ^-^e the modified Bessel functions of the first and second kind. The 
term Hd is the solution of the homogeneous London equation 

Hrf + rotrotH^ = (7) 

in the half-space z <0 with the boundary condition d). To solve equation (|^) it is conve- 
nient to represent the components H^^ = {H^, H^, H^) in the form of the two dimensional 
Fourier transform over spatial harmonics: 



^ +00 +00 

H2 = — j dqexp^iqx) J du exp (iuy) (q, u) exp {z\Jl + + v?), (8) 

— oo — oo 

a = X, z. 



Equation (0) then reduces to the simple algebraic relation for the coefficients C2{q,u): 



-qC^^-uC',+z^l + q' + u'C^, = 0. (9) 

In compliance with (H), the magnetic field can be written as a negative gradient of a 
scalar potential f/^: 

H, = -WUs, (10) 
and the potential Us obeys Laplace's equation 
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AUs = 0. 



(11) 



It is evident, that the solution of equation ([Tl|) can be represented in the form 



^ +00 +00 
= — J dqexp{iqx) J du exp {iuy)Cs{q, u) exp {—z\J + u^). (12) 

—00 —00 

The boundary conditions for the magnetic field (0) lead to three algebraic equations for 
the so far undetermined coefficients C^(g, m), Cs{q,u): 



Cl = -iqCs, Cl = -iuCs, Cl--SM,u) 



g2 + m2 



(13) 



The solutions of the linear algebraic system (|9|, [T3|) for + 7^ are 



qJl + q^ + S^{q,u) 



nJ q'^ + u'^ I Jl + q"^ + u"^ + Jq"^ + u'^ 



uJl + q^ + v? Sv{q,u) 



nJ q^ + \ Jl + q^ + + J q^ + 



i\Jq^ + Sy{q,u) 



n[^l + q^ + u^ + \/q^ + u^ 



+ q^ + Sy{q,u) 



irJ q'^ + u'^ ( a / 1 + g2 + + -i/g^ + 



(14) 



If + = 0, then the coefficients C^, Cg in the expansions (||,0) are equal to zero. 
The function Sv{q,u) is determined by the spatial Fourier spectrum of the magnetic field 
distribution Hy in the plane z = 0: 



, -f-OO -f-OO 

Sy{q,u) = J dxexp{-iqx) J dy exp {-iuy)Hy{y, x) 



(15) 



Substitution of (|) in the expansion (0) leads to the following representation of the 
spectral function Sy{q,u): 



N • / ^ ) hiyy^ + (i^)Kiipv\^'^ + Q^)^ y<pv; 

Sy{q,u)=py / dysm{uy)< " " ^ 

i Ui(p,v^rT7)Ki(i/v^rT7), y>p,. 



(16) 



The details of the numerical simulations of the spectral function Sy{q,u) are presented 
in Appendix A. Using the asymptotic expansion of the modified Bessel functions /i(C); 
Ki{() for C ^ 1 one can obtain 



Sy{q,u) 



sin (upy 



u 
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M > 1, q, 



(17) 



which is a good approximation of the exact expression (p!6| ) at the large values of spectral 
variable u. In the other limiting case ( g,M <^ 1, wp^, ^ 1 ) it is convenient to represent 
the spectral function Sy{q,u) (p!6|) as 



Sy{q,u) c::^ au, g,M-Cl, wpt, <^ 1, 



where the factor a is expressed in terms of the vortex radius p„ : 



a = Pv 



Pv +0O 

Kr{pv) j dCCh{0 + h{Pv) I dCCK^{C) 



Pv 



(19) 



Thus, the magnetic field structure of a solitary vortex half-ring placed near the flat 
surface of a type-II superconductor is fully determined by the relations (^J^J|,|^,|^,|l^,|l^) . 
The field describes the disturbances of the structure of the toroidal Abrikosov vortex 
(H) by the surface S. In the area of the vortex core the field is in opposite direction 
to the vortex line, so that (H^ ■ H^,) < 0. It results in a decrease of the magnetic field 
value at the VL. The superconductor surface effect is strong only in the surface layer of 
a thickness defined by the vortex radius p^ or the magnetic field penetration depth A. 
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The field denotes the scattered magnetic field outside the superconductor, pro- 
duced by the isolated fiux line in the form of the half-ring placed near the surface. To 
determine a structure of the field let us estimate a behavior of the potential Us (0) 
at long distances 2; ^ 1, p^, from the surface. Since the integrand in (|l^) drops abruptly 
if q"^ +u'^ > (1/z)^, only small values of the spectral variables q and u are important. As- 
suming that q, u <ti I, 1/ Pv and using the approximate expression of the spectral function 
Sy{q,u) (0), the potential Us ([l^) can be written in the form: 

TT ^0. f 1 / \ f 7 usm[uy) . rz , , 

~ — / aqcos{qx) / au . ^ exp (— 2:^/^"' + w ), z^l,p.o. (20) 

yQ -ru 

It is shown in Appendix B, that the potential (|20|) is the potential of the dipole P = 2p^Q 
(see the inset of fig. H), if the effective charge Q = a/27rpt,. Figure plots the values of 
the charge Q and dipole moment P versus the vortex radius py. The dipole structure of 
the potential Us determines the asymptotic decay law for the scattered field at long 
distances z ^ 1, p^ from the surface: Hs ~ R~^, where R = {x'^ + y"^ + z^'f^'^ ^ p„ is the 
distance from the center of the vortex half-loop. 

To calculate the magnetic fiux $ in an Abrikosov vortex, the cross-section S has 
to be selected perpendicular to the vortex line. In the case of the vortex half-ring the 
intersection plane E has to include the center of the semicircle formed by the VL and 
the axis x (see fig. |l|). The specific feature of magnetic structures in the form of a vortex 
half-ring is self-consistent localization of both the magnetic field and supercurrent in the 
vicinity of the VL. Since the magnetic field and the supercurrent vanish at large distances 
from the VL, the total magnetic fiux $s across the plane S is equal exactly to zero: 

<l>s = $E_ +$E^ =0, (21) 

Here $s_|_ are the magnetic fiuxes across the corresponding half-planes Ej-. Using the 
Gauss theorem, it is easy to verify that the fiuxes do not depend on the orientation 
of the intersection plane E relative to the surface of the superconductor. So, by choosing 
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for simplicity the plane y = as the cross-section S, one can obtain: 



(22) 



-oo 



$^ = J dx J dzH^{z,x) 



+00 

$rf = dx dzH^{x,0,z), 



-oo — CXD 



-OO — OO 



hoo -l-oo 



$2^=$^= y rfx y dzHy{x,0,z) 



(23) 



-oo 



The term is the magnetic flux in the solitary toroidal Abrikosov vortex and was 



calculated in [11 



(24) 



and the magnetic fluxes $d and $s are determined completely by the spectral function 



$d = / du 7===, 



(25) 



^^^-'-Tdu^^l^^. (26) 
J u[u + V1 + ?? 





Using the relations (^,^,^) and the definition of the spectral function Sy{q,u) (0) it 
is easy to verify that the total magnetic flux $s is really equal to zero. The results of the 
simulations of the magnetic fluxes $i, and $s_ versus the vortex radius p„ are shown in 
flg. ^ The magnetic flux $s_ of the vortex half-ring depends strongly on the VL radius 
Pv and tends asymptotically to the quantum of magnetic flux $o at large values of p^. The 
flux $s_ differs noticeably from the flux of a toroidal Abrikosov vortex $i, . It means 
that the surface of a superconductor affects signiflcantly the magnetic fleld structure of 
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the vortex half-ring, thus largely reducing the value of the magnetic field in the vortex 
core area [0. 

In the region z < 0, occupied by the superconductor, the London equation (|^) 
implies the following expression for the free energy G- 



G_ = — 
Svr 



J rfx c/y [h X rot h] + J rfVH(H + rotrotH). 



(27) 



z = 



z<0 



The energy contribution of the magnetic field outside the superconductor {z > 0) is 
defined in the usual way: 



z>0 



{2i 



The total free energy Go is the sum of the expressions ( pTlpH]) : 



Gq — G— -|- G_| 



(29) 



To find the free energy of the vortex half-ring it is necessary to substitute the solutions 
(^|6|,|8|, [T0|JT^ , P^JT6[) describing the magnetic field distributions H and in the expressions 



( P7| , PH| , PUD . It is convenient to extract the contribution of the superconductor surface and 
represent the energy Gq as the sum: 



Go — + Gs- 



(30) 



Here the first term G^ is the free energy of the half of the toroidal Abrikosov vortex. The 
magnitude of G^ is fully determined by the value of the magnetic field in the center of 



the vortex line [10, 



Pv 



Gv = -—Hy{pv, 0). 



(31) 
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The second term Gg includes the total influence of the surface S and is expressed in terms 
of the spectral function S^i^q, u): 



G 



hoo +00 



2tx 



3 



dq / du ' . , X , (32) 





7r/2 



D(q,u) = — — ^ — [ dif cos If cos (upv simp) 
up,, J L 







2 I 2 
q + u 



sin (y9 sin {up^ sin (p) exp {—pv cosifJl + + u'^). 



U\/\ + q^ + u^ 



Expression (|^) has a logarithmic divergence at the VL center (p = p^,, x = 0), since the 
London equation (^) is not valid in the region of the normal vortex core. To compute the 
energy G^ (pi]), the magnetic field H^{py,0) was thus truncated at the distance C, from 



the VL center. The details of the numerical simulations of the magnetic field H^{p^,0) 
are presented in Appendix C. Figure ^ gives the dependencies of the energies G^ and Gs 
on the radius of the vortex half-ring p^. Let us note that Gs ^ at any values of the 
vortex radius. It means that the free energy Gq is modified insignificantly by the surface 
of the superconductor and the estimate 

Go(p.„) ^ GM (33) 

accurately approximates the exact value of the free energy of the vortex half-ring. 

Thus, the free energy Go(Pd) of the Abrikosov vortex in the form of a semicircle 
of radius p^ exceeds very slightly the energy G^{pv) of half of the toroidal vortex. At 
the same time, the magnetic field structures of these vortices differ noticeably, especially 
inside the layer with thickness ~ A near the surface. The surface changes the structure 
of the toroidal vortex, such that the value of the magnetic field in the core of the vortex 
decreases. As a consequence, the magnetic fluxes (11)23) of the vortex half-ring 

are smaller than the magnetic flux (^4]) of the toroidal Abrikosov vortex even if the 
size of the vortices exceeds the magnetic field penetration depth A (p„ > 1). Due to the 
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continuity of the magnetic flux, the flux line creates the scattering field outside the 
superconductor. The energy of the magnetic field H,, makes an important contribution 
to the total free energy of the vortex half-ring and neutralizes the decrease of the London 
energy (13) > forced by the distortions of the vortex structure, induced by the surface 
of superconductor. 



3. VORTEX HALF-RING IN AN APPLIED MAGNETIC FIELD 

The behavior of a superconductor in an applied magnetic field Hq is controlled by 
the Gibbs free energy [|]: 



G-^ I dV{HHo). (34) 



47r 

2<0 

If the applied field Hq is directed along the surface S, the magnetic field distribution in 
the superconductor half-space {z < 0) (see fig. |lD can be represented as follows: 

H = Hy + H^. (35) 

The first term Hy obeys equation (|l|) with the boundary conditions (^ and describes 
the magnetic field structure of the Abrikosov vortex with arbitrary vortex line, if the 
applied field is absent. The second term = iioexp{z) is the Meissner (vortex- 
free) distribution of the magnetic field near the surface of the superconductor. The free 
energy G is defined by the expressions (^,^,^), where the magnetic field H in the 



superconductor satisfies relation (|35|). Taking into account that the field on the 
superconductor surface S coincides with the applied magnetic field Hq, the Gibbs free 
energy F (0) measured relative to the free energy of the Meissner state can be expressed 
as: 



F = Gv + ^ y ^^yiim ~ ^ J ^^vHo 



(36) 
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Here Gy is the free energy of an Abrikosov vortex described by the term Hy in (^) in 
a zero apphed magnetic field. The VL location of this vortex in space is determined by 
arbitrary curve 1^ and the vector = dl^By is directed along the vortex line. It should 
be noted that expression (|36D sets the relation between the free energy of the Abrikosov 
vortex Gy and the Gibbs free energy F for arbitrary shape of the vortex line. Using (|36|), 
it is easy to write the Gibbs free energy of the Abrikosov vortex in the form of a semicircle 
of arbitrary radius p^: 

tt/2 

F{py) = Go{py) - 1 - J dipsmipexp{-p^smip) . (37) 



Figure |^ gives the dependencies of the Gibbs free energy F (^Tj) on the radius of the vortex 
half-ring p^ for some values of an applied magnetic field Hq. The type of the curve F{py) 
depends strongly on the value of the applied magnetic field Hq. Taking into account the 



equality (33) and considering that in the case p„ ^ 1 the following approximate expression 



for the free energy G^{p^) holds 



Gv{Pv) - —^Pv, Pv > 1, (3^ 



one can obtain the asymptotic behavior of the Gibbs free energy F{p^) for p^ 3> 1: 

p. > 1. (39) 



— Hp^ — Hn 

2 



Here Hd = In {k)/A7t is the lower critical magnetic field for the London model. If Hq < 
H*^ = 'kHci/2, the Gibbs free energy F{py) is positive for any p^ ^ 0, and the Meissner 
state has the minimal energy. So, the existence of vortex half-rings in the superconductor 
is forbidden, and only vortex-free distributions of the magnetic field and supercurrent are 
possible. In the case of Hq = H*^, the Gibbs free energy of a state containing a vortex 
half-ring of a large radius p^, ^ 1 matches the energy of the Meissner state. It should be 
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noted that the vortex half-ring is asymptotically stable in this case, because 

dF 



lim 



0. 



So, the field H*^ is the minimal applied magnetic field, when it is energetically favorable 
for a vortex half-ring to be formed near the fiat surface of type-II superconductor. If 
Hq > H*^, there is a critical radius pc of the vortex half-ring at which the curve F{py) has 
its maximum U (see the inset in fig. Figure ^ gives the dependencies of the height of 
the surface energy barrier U and the critical radius pc on the value of an applied magnetic 
field Hq. Since the Gibbs free energy F^p^) decreases monotonically at p^ > pc, the 
energy barrier U prevents nucleation of expanding vortex loops near the surface. The 
barrier height U is reduced by increasing the applied field Hq. When the applied field Hq 
is large enough for the critical radius pc to become of the order of the coherence length ^, 
the surface barrier U vanishes. Using the approximate expression for the free energy G„ 



of the vortex half-ring |jT5|, [16 1 



In(Pf) 
Air 



pv <l 



(40) 



one can estimate that the surface barrier U is destroyed 

dF^J__pJlQ 

dp, - IGtt 4 ^' 
when the applied magnetic field Hq becomes of the order of the thermodynamic critical 



field in the superconductor He 



K 



271 V2 



. Thus, the same applied magnetic fields are 



required to destroy both the surface barrier for a vortex half-ring and the Bean-Livingston 
barrier for a rectilinear VL parallel to the surface. 



4. DISCUSSION OF THE RESULTS 

The solutions of the London equations were obtained here to describe the structure 
of a curved Abrikosov vortex near the fiat surface of a type-II superconductor. The VL of 
such a vortex is a semicircle of arbitrary radius. The interaction of this vortex half-ring 
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and an applied magnetic field was considered by calculating the Gibbs free energy. This 
problem extends the well-known ideas concerning the surface Bean-Livingston barrier [^], 
for the case where the Abrikosov vortex is curved. Vortex loops like these represent an 
example of compact magnetic structures in SC fl^. The specific feature of such magnetic 



structures is the self-consistent localization of both the supercurrent and the magnetic 
field in all three dimensions. Due to the continuity of the magnetic flux, the end points of 
the vortex loop touch the surface of SC, and the vortex line produces a scattered field 
outside the superconductor. The creation of the vortex loops is depressed significantly by 
the surface energy barrier. The height of this barrier is reduced by increasing of the applied 
magnetic field o- The surface barrier disappears completely when the field o becomes of the 
order of the thermodynamic critical field Hem-, and the deterministic nucleation of small 
vortex loops with size comparable to the coherence length ^ is possible. The applied 
magnetic field Hq produces the Meissner supercurrent 



-^rotH^, (41) 
47r 



distributed mainly in the layer of thickness A near the surface. At the surface the density 

of the Meissner current, corresponding to the applied field Hq = Hem, is comparable to 

the depairing current density = ^ „ 0. Since the inhomogeneity of the Meissner 

12v 3 TT 

current density is not important at scales of the coherence length ^, the condition jm ~ jc 
determines the threshold of thermodynamic stability of the uniform current carrying state 
in superconductors ^ and agrees with the Landau criterion of roton excitation in a 
superfiuid flow 1T|, 0. 

At a nonzero temperature T, before the surface barrier is destroyed {Hq < Hem.), the 
penetration of magnetic flux into the SC occurs due to thermally activated nucleation, 
subsequent growth, and merging of vortex half-rings. A similar mechanism of thermally 
assisted overcoming of energy barriers should be important at high temperatures T near 
Tc [0 especially in materials with strong thermal fluctuations [|1^, pO[] . 



The creation of expanding vortex loops near a superconductor surface is similar to 
the generation of closed Abrikosov vortices by an external current |l^, |T^. Since the 
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energy of the vortex half-ring Gy (31) grows monotonically with the radius of the VL 
p^, the nucleation and expansion of vortex loops results in the increase of the total free 
energy of superconductor. The decrease of the Gibbs free energy F (0) is determined by 
the work AH^(p„), performed by the source of the applied magnetic field Hq during the 
expansion of the vortex loop up to size p^: 

■k/2 

AW{pv) = 1— / d(p simp exp {— py sin (f) . (42) 

27r L 7 J 



The work AW{py) is related to the motion of the vortex line under the action of the 
Lorentz force /l 

h = [jm X e„] (43) 

caused by the Meissner current jm- Since the current decreases outside a surface layer 
of thickness ~ A, the interaction of the applied magnetic field and the vortex half-ring 
depends strongly on the size of the vortex py. If p„ <^ 1, the vortex line is located wholly in 
the area where the current density is distributed uniformly, and AW ~ p^,^ |T^, |T^]. 
In the other limiting case (p^, ^1), the integrand in the expression (|4^) drops abruptly if 
sin V? ~ 1 , and the region sin (p < l/p^ <^ 1 makes an essential contribution to the integral 
in (53). Replacing sin 99 ~ 99 in (H2|) it is easy to obtain 



AW{py)c^^ for p,»l. 

ZTT 

In this case the basic part of the vortex line is located outside the surface layer of a 
thickness ~ A, where the Meissner current jm is absent, and the work 6W is produced by 
the motion of parts of the VL, that border on the surface. 

The equilibrium shape of the Abrikosov vortex is determined mainly by the structure 
of the applied magnetic field and represents a rectilinear VL parallel to the surface. The 
semicircle is the correct form of a vortex nucleus near the surface of the superconductor 
only if the vortex radius p^, ^ 1. In this case the Meissner current jm does not change 
essentially on scales of the size of the vortex, so the value of the Lorentz force (|^) is 
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the same for different points of the vortex hne. Thus an expansion of the vortex half-ring 
takes place without distortions of the form of the vortex line. When the size of the vortex 
Pv becomes comparable to the magnetic field penetration depth A (p„ > 1), the value of 
the Lorentz force fi reduces strongly for the parts of the vortex line, placed farther from 
the surface. It results in distortions of the form of the vortex line so that vortex loops 
elongate parallel to the surface. It is evident that the vortex line of such an elongated 
vortex loop has a large part directed along the applied magnetic field. 

The pattern of expanding single vortex half-rings holds as long as the radius of 
the vortex is much smaller than the typical intervortex distance L. If the intervortex 
distance becomes comparable with the sizes of the vortex, the regular expansion of vortex 
half-ring is changed by the interaction of adjacent vortices. The attraction of the sections 
of the vortices of opposite directions may lead to intersection and cross joining of the 
adjacent vortices 0. Figure |^ shows the sequence of formation of a rectilinear Abrikosov 
vortex parallel to the surface by joining of expanding vortex loops. It is essential, that if 
the applied magnetic field is not enough to destroy the surface barrier {Hq < Hem), and 
the thermally activated motion of vortices over the barrier takes place, the creation of 
a rectilinear Abrikosov vortex near the surface by intersection and joining of the vortex 
loops can allow for the magnetic field entry in type-II superconductor. The surface barrier 
in this case is determined by the energy barrier U preventing generation of a single vortex 
half-ring and may be lower than the ordinary BL barrier for a rectilinear Abrikosov 
vortex IQ. The adjacent vortex loops, which penetrate after the first one, require less 
energy than U fl^. Therefore, the process of the formation of a rectilinear Abrikosov 



vortex is of avalanche-type, and U is the activation barrier for the vortex entry . 



5. CONCLUSION 

In this paper the properties of the Abrikosov vortex in the form of a semicircle of 
arbitrary radius p„ near the fiat surface of type-II superconductor in the applied magnetic 
field Ho, parallel to the surface, have been theoretically studied. The magnetic field 
structure, the magnetic flux and the free energy of a vortex half-ring have been calculated 
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in the London approximation. The creation of such vortices is delayed by the surface 
barrier, similar to the Bean-Livingston barrier for a rectilinear Abrikosov vortex [|]. 
The barrier is reduced by increasing the applied magnetic field Hq and vanishes in the 
thermodynamic critical field Hem, so that the deterministic nucleation of vortex loops 
with the size of coherence length ^ becomes possible. At temperatures T near and 
for applied fields Hq < Hem the magnetic flux penetrates into a type-II superconductor 
by thermally activated nucleation of expanding vortex half-rings near the surface. The 
formation of the rectilinear Abrikosov vortices by intersection and cross joining of vortex 
loops has been considered. 
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APPENDIX A 

Let us express the spectral function Sy{q,u) by means of the auxihary function 
Wqu{z), satisfying the following inhomogeneous differential equation with the zero initial 
conditions: 



1 , 



- -'^qu - (1 + 9 )wgu = - sin (uz), {A.l] 



Writing the Green function for equation (A.l), one can obtain the particular solution: 



Wqu[Z) 



-- zK^{z^Jl + q^) J dC sin «) I,{C^l + q^) 



2 

-zh{z^l + q^) J dC sin«)Xi(cVl + g')- 



Using the known integral: 

+00 

/7iU 
dCKi{aC) sin«) = 

one can set the following relation between the functions Sv{q,u) and Wgu{z): 

Sv[q, u) = WguiPv) + I / (A.2) 

2V 1 + ? V 1 + ^ + ^ 

Thus, to simulate the value of the function S'„(g, u) for arbitrary q and m it is possible to 
solve the equation (A.l) on the interval z G [0,Pt,] and to substitute the obtained value 
Wqu{py) in expression (A.2). This method of calculation of the spectral function Sy{q,u) 
is especially efficient for large values q,u^ 1, because it uses the well-known algorithms 
of solution of Cauchy problem instead of simulations of improper integrals on rapidly 
oscillating functions. 
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APPENDIX B 



The potential Up of the dipole P = 2p^Q placed in the plane z = as shown in the 
insert of fig. ^ is described by the well known expression 



Up = P^, R={x' + y' + z'f':>p,. (B.l) 

Let us represent the potential Up in the form of the two dimensional Fourier transform 
over spatial harmonics: 

^ +00 +00 

Up = — / dqexp{iqx) / duexp {iuy)Sp{q,u, z). (B.2) 
zvr J J 

—00 —CO 

The spatial spectrum Sp obeys the inverse Fourier transform and should be written as 
follows: 

^ +00 +00 
Sp = — / dxexp{—iqx) / dy exp {—iuy)Up{x,y, z). (B.3) 

—00 —00 

Substitution of Up (B.l) in the expansion (B.3) leads to the following expression of the 
spectrum Sp: 

2iP Y I 

Sp = J d( cos K) Ko{q ^Jz^ + e), (B.4) 



where Kq{() is the McDonald function. Using the known integral: 



J dC cos{uOKo{q^z' + e) 




TT exp {—zJ + u'^) 



2 Jg^ + 



one can obtain: 



u exp {—zJ q^ + u^'^ 



Sp{q, u, z) = -iP ^ . (B.5) 

Thus, the potential Up of the dipole P at long distances R ^ is determined by the 
expression: 

IP f , \ [ 1 usmiuy) rz 

Up = — / dqcos{qx) / du . ^= exp {—zJ q^ + w'). (fj.o) 
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APPENDIX C 



To compute the magnetic field value at the vortex line center, we use the other 
representation of the magnetic field distribution Hy{p,x) in a toroidal Abrikosov vortex 



IC 



where Ji(C) is the Bessel function of the first kind. It is easy to verify, that expression 
(C.l) is the exact solution of equation (^. The integral in (C.l) is divergent for p = py, 
X = 0, because of the integrand function falls down too slowly. To estimate the magnetic 
field value Hy{py,0) let us extract the divergent part of integrand function 

Sh{q,p) = 7rT=^ • 

Using the asymptotic expansion of the Bessel function Ji(C) for ^ 1, one can obtain 
the following behavior of the function Sh{q,p)'- 

The first term in (C.2) results in the logarithmic divergence of expression (C.l) if p — > p„, 
a; — > 0, because of 

cos[g(p-p^)] 

/ /r^-^ =^o{\p-Pv\), C.3 

where Kq{() is the McDonald function. The value of the last integral have to be truncated 
in the usual way at the distance ^ from the VL center by setting |p — p^| = ^ in (C.3), so 
that the following representation of the magnetic field value if^(p^,0) at the VL center 
holds: 



u ( n\ ^oiO - 7r/2 p„ . 
Hvipv, 0) = + — dq 



qJl{,qPv) 1 ^ sin(2gp^) 

y/lT~q^ vrp^Vl + '^QPv 



(C.4) 



The new integrand function in (C.4) falls down quickly enough to provide the convergence 
of the improper integral in expression (C.4). 
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Figure 1: Abrikosov vortex in the form of the semicircle of the radius p,^ near the flat 
surface S of type-II superconductor. The superconductor occupies the half-space ^ < 
and is characterized by the magnetic field penetration depth A and coherence length ^ ; the 
vortex line is placed in the plane {yz) . The intersection of vortex plane and cross-section 
E = E_ + E+ is shown. 
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Figure 2: Effective charge Q (1) and dipole moment P — 2pyQ (2) vs. radius of the vortex 
hne py. The inset shows the dipole P placed in the plane z — 0. 
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Figure 4: Free energy of half of toroidal Abrikosov vortex (1) and tiie contribution Gg 
of the surface of superconductor S to the free energy of the vortex half-ring (2) vs. the 
radius of the vortex hne {k — 100). To maintain the scale of the vertical axis fixed, 
the value 100 • Gg is used . 
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Figure 5: Gibbs free energy F of the vortex half-ring vs. the radius of the vortex hne 
p„ for several values of an applied magnetic field Hq {k — 100). The numbers near the 
curves denote the value of the field Hq in the units of the lower critical field Hd- 
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Figure 6: Surface barrier U (1) and critical radius of tiie vortex iialf-ring pc (2) vs. the 
value of the applied magnetic field Hq [k = 100). The inset sketches the part of the 
curve F{p^). The Gibbs free energy F{py) has its maximum U at the critical radius Pc- 
F{pc) = U. 
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Figure 7: The formation of the rectihnear Abrikosov vortex, parallel to the surface, by 
intersection and cross joining vortex loops. 
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